Development of microsatellites markers for the deep coral Madracis myriaster (Pocilloporidae: Anthozoa)

In 2013 Colombia made an important step towards the construction and management of Marine Protected Areas (MPAs) by establishing the first Deep Corals National Park (PNNCP). Inside this MPA, the coral Madracis myriaster (Cnidaria: Pocilloporidae) was found as the main reef builder, offering habitat for many species of fish and invertebrates. In order to improve the study of deep-sea coral habitats, their connectivity and prospective management, nine new genetic markers (microsatellites) were developed for M. myriaster and tested in samples from PNNCP. We present the assessment of these markers, with a specificity for the deep coral, and its prospective use in future analysis for the PNNCP and other areas in the Caribbean and the Atlantic, where M. myriaster is reported. We also include an additional taxonomic analysis performed on samples of M. myriaster using scanning electron microscopy.

PCR standardization for microsatellites and DNA sequencing. The test of the 14 microsatellite loci obtained from the filtering process (see methods) with the 26 M. myriaster samples showed correct amplification in 13 samples. The other 13 samples had clear electrophoresis bands, but unclear reading of the fragment sizes in "Fragman" and poor binning in 'MsatAllele' (Table 1). A total of nine microsatellite loci were confirmed for the deep-sea coral M. myriaster, using the analysis in "Poppr" and the criteria of 20% of missing data across the 13 samples. The primer specificity test, using samples from M. auretenra, Montipora sp. and Antillogorgia sp. under the same amplification conditions used for M. myriaster, confirmed no amplification on the other species. Data analysis. The assessment of the genotyping results, performed with the final 26 DNA selected samples, showed the potential of these loci for further population genetics analyses. The test for clones, based on multilocus genotypes ('Poppr'), did not find clones in the samples used. The allele number per locus (Na) was between 4 and 11. The observed and expected heterozygosity ranged from 0.334 to 0.81 and 0.66 to 0.88, respectively

Discussion
This study presents nine new microsatellites developed for M. myriaster, based on the samples collected at PNNCP (Colombia). All the new microsatellites were tested in all the available samples, independently of their quality of DNA. Despite obtaining microsatellite amplification of 14 loci in the 26 samples of M. myriaster, with clear bands in the electrophoresis agarose get prior to the genotyping, 13 samples were excluded because some of their fragment sizes under "Fragman" software did not reach the threshold to detect them, or because the amplification was not clear to bin the fragment sizes using "MsatAllele". Understanding that the exclusion of www.nature.com/scientificreports/ that data limited the study scope, the decision was taken to keep a conservative approach and eliminate false positives in the analysis. Likewise, five microsatellite loci with low amplification in the 13 samples were excluded from the analyses. These decisions were taken to avoid possible data misinterpretation from amplification errors caused commonly by low concentration of DNA from samples kept in museums 31 , as was the case of these 13 samples with a DNA concentration below of 10 ng/μl. The low DNA concentration could create a common error during the PCR amplification, called "false homozygotes", due to random sampling in the DNA template and deficiency amplification as a response to very low DNA quantity 32 . The literature suggests increasing the number of PCRs to reduce the occurrence of "false homozygotes" 32 ; unfortunately, it was not possible to perform the additional PCRs due to lack of additional DNA material from the samples. However, it is highly likely that any of the 14 developed microsatellites, including the last five removed markers, can be amplified using samples with better quality and quantity of DNA.
The newly developed microsatellites for M. myriaster may be used as a tool to examine genetic variation among populations of M. myriaster in the Caribbean Sea once more samples from the area are available. The usefulness of these type of molecular markers in conservation studies has been reported 19 and they have been used successfully with shallow corals 33,34 and deep-sea corals 27,28 . In addition, microsatellites have also proven to be helpful in the design and effectiveness of MPAs, based in the data given for populations of connectivity and genetic structure derived from these markers 35 . Molecular markers have also been incorporated, as complementary tools, in the identification of species in the genus Madracis. Some examples include the Internal Transcribed Spacers (ITS) for genetic variation 13 , the putative mitochondrial control region-CR 36 , the ATP8 mitochondrial marker 37 , as well as the ATPS and SRP54 nuclear markers 38 . For this reason, the specificity of the new microsatellites developed for M. myriaster shows the potential use of these markers as a complementary taxonomic identification tool, because it was found that they amplify only in samples of this particular species and no other corals.
In this study, nine samples catalogued as M. cf. myriaster were excluded at the microsatellite standardization stage due to the non-amplification, despite the good quality and quantity of DNA. These samples could have been M. brueggemanni due to the morphologic similarity with M. myriaster 2 . However, it was not possible to confirm the identity of the samples catalogued as M. cf. myriaster, because they had the same morphological features under traditional taxonomic methods but were identified as "different" by the microsatellite standardization method. Also, the rudimentary method used to collect the samples in most of the stations in Colombia, known as trawl net sampling, causes difficulty for the correct identification of the fragments because the method breaks and mixes the colonies from different species into small fragments over a 1 km transect sampling 12,39 . These small coral fragments (less than 2 cm) are very difficult to identify under traditional taxonomic methods, as it was evident in the results, because there are not significant differences in the branch thickness, size or, the corallite diameter among samples, characters commonly used in traditional taxonomy.
Despite the obvious limitations to perform additional analyses due to the reduced number of samples currently available, additional constraints regarding the high costs and logistics to obtain deep-sea coral samples obstruct the exploration on these habitats. However, the new microsatellites presented here can be used in future analysis, considering the need to continue the exploration of marine national parks, MPAs in Colombia, and other areas in the Caribbean. Currently, the PNNCP has only been explored until 500 m depth, but there is a need to explore the deepest areas of the park, reaching 1350 m depth 40 . Personal communications with other colleagues suggest that this coral could also be present at similar depths in other Caribbean and Atlantic locations (e.g. Gulf of Mexico, Florida and Brazil) 3,41,42 . Therefore, these markers may be used in the near future as a tool to strengthen the incomplete identification of samples, to complete the base line of deep environments in MPAs, to examine the Caribbean connectivity among deep-sea habitats, and to develop comprehensive conservation management plans of MPAs in the Caribbean and Atlantic region.

Methods
Sampling. Fifty-three fragments of M. cf. myriaster were used for this study from the reference collection of Cnidarians at the Marine Natural History Museum of Colombia (MHNMC), part of INVEMAR. The samples belong to different exploration projects in the Colombian Caribbean Sea, from 1998 to 2015, where trawl net methodology was used at depths between 120 and 350 m, following transects of 1 km in length. After collection, all samples from the same transect were mixed and preserved in 70-90% ethanol. Table 2 shows the initial coordinates of each transect, from where samples were selected for molecular analysis (more information: https:// siam. invem ar. org. co/ campa nas-proye ctos). The 53 samples were exported to Manchester Metropolitan University under the CITES permits 40,885 (25th July 2016), 41,449 (25th January2017) and 43,908 (6th May 2019).
Analysis on Madracis myriaster samples. All the 53 coral fragments were originally examined to confirm the correct identification because the method used to collect the samples in Colombia, known as trawl net sampling, usually breaks the colonies into small fragments, mixing the fragments from different colonies over the 1 km transect sampling, and increasing the chance of mixing colonies of M. myriaster and M. brueggemanni. However, 18 samples originally catalogued as M. myriaster were eliminated due to degraded DNA (DNA extraction, PCR standardization and DNA sequencing section). Another nine samples originally catalogued as M. cf. myriaster were eliminated from the analyses due to unsuccessful amplification (DNA extraction, PCR standardization and DNA sequencing section). As a result, the identification effort was focused on the 26  www.nature.com/scientificreports/ myriaster or "fluted finger coral" often has branched colonies with an irregular morphology and purple, pink or orange colours. In contrast, M. brueggemanni forms colonies with tiny and delicate branches, giving a less complex three-dimensional structure arranged in three dimensions. M. brueggemanni also presents a narrow distribution across the Western Tropical Atlantic and is found between 51 and 160 m deep 2 . Twenty-four from the 26 samples were further examined due to ambiguous morphology using scanning electron microscopy-SEM (Zeiss Supra 40VP FE-SEM with an EDX detector for elemental analysis at Manchester Metropolitan University -MMU), due to the mismatch with the species descriptions above and the extremely small size of the fragment for the comparison. The images were analysed using the program Image J to quantify branch thickness (three measurements per photo), corallite diameter (two measurements per photo) and presence of spines in the sample. Afterwards, the SEM images from the 24 samples were revised with Dr Nadia Santodomingo, at the Natural History Museum of London-NHML, splitting the samples in two groups (group a with 11 samples confirmed as M. myriaster and, group b with 13 samples as M. cf. myriaster) in order to find morphological differences using a t-test.
DNA extraction, PCR standardization and DNA sequencing. The following DNA extraction protocol was optimized using two samples identified as M. myriaster (samples codes 1619 and 1620) with the DNA tissue and blood extraction kit from Qiagen: (1) fragments of the sample (0.5-1 cm diameter) were dried at 36 °C for 10 min using a Thermomixer; (2) 180 μL of buffer AL and 20 μL of Proteinase K were added to the sample and vortexed for 20 s; (3) the samples were incubated in the Thermomixer at 55 °C for 2 h with vortex every 30 min; (4) after digestion, the fragment was removed (calcium carbonate skeleton) and 200 μL of buffer AL was added, after a brief vortex for 20 s, 200 μL cold ethanol was added; (5) all the mix was transferred to filter tubes and centrifuged at 10,000 rpm for 5 min; (6) the filter from the tube was transferred to a new collection tube, 500 μL AW1 were added and centrifuged at 10,000 rpm for 3 min; (7) the filter was transferred after the centrifuge into a new collection tube, 500 μL AW2 were added and centrifuged at 13,500 rpm for 3 min; (8) the filter tube contents were transferred into a 2 ml Eppendorf tube and 20 μL of buffer AE (pre-heated at 75 °C) were added; (9) then left for 15 min at room temperature (10) afterward centrifuged at 10,000 rpm for 1 min; (11) steps 9 and 10 were repeated and then the filter was discharged; (12) the final sample (2 ml Eppendorf tube) was stored at -20 °C. Electrophoresis was used to test the quality of the DNA extractions (53 samples) with the following conditions: (1) 0.80 g of agar were mixed with 80 ml of 1X TBE buffer solution (1% agar gels) and 0.8 μL of GEL GREEN to prepare the gels; (2), 1 μL of sample plus 2 μL of loading buffer were added into each pool. The electrophoresis chamber was run at 60 V for 60 min. At this point, 18 samples confirmed as M. myriaster were excluded due to absence of DNA after the extraction. The DNA from the two samples chosen to improve the extraction protocol (1619 and 1620-good quality and quantity of DNA) were normalized to 50 ng using the Nextera® DNA sample Preparation Kit and sent to the sequencing facility at the University of Manchester-UoM (UK) to perform an Illumina MiSeq paired-end sequencing (2 × 250 bp). The sequencing data were analysed using the bioinformatics tool Palfinder on the Galaxy Centaurus Server Platform (https:// palfi nder. ls. manch ester. ac. uk), and the primers were designed using the optional filters pipeline, with parameters for melting temperature, annealing temperature and primer length of the Type-it Microsatellite PCR QIAGEN kits 18 . From the process, 223 potentially amplifiable loci (PALs) were obtained and 36 PALs (7 tri-and 29 tetra-nucleotide motifs) with good motifs and GC contents were selected for further analysis. Nine samples identified as M. cf. myriaster were eliminated at this point due to unsuccessful amplification. The genotyping was performed with the remaining 26 samples ( Table 2). The fluorescence labelled PCR products were sent to the University of Manchester DNA Sequencing Facility (UK) and to the Core Genomics Facility at the University of Sheffield (UK), in a mix of: 9 μl of HiDi Formamide, 0,2 μl of the Liz 500 (GeneScan™ 500 LIZ®) and 0,8 μl of the PCR product (including the positive and negative controls). The products were sized in both places using the capillary electrophoresis Applied Biosystems 3730 DNA Analyser (enabling size discrimination within the range of 20 to 600 base pairs using a range of dyes). A species-specificity primer amplification test was also performed on M. auretenra ( www.nature.com/scientificreports/ genotyping. The values of each allele size were normalized using the positive control. The 'MsatAllele' package 47 was employed to bin the fragment sizes (Fragman output), assigning to a set of defined alleles the closest size for each marker. Then the true allele set was determined using the known repeat length and using the histograms of observed fragment lengths before and after the binning. The package 'Poppr' in R 48 was also used to assess the frequency of missing alleles across primers. Additionally, the number of clones present in the sampling was investigated by counting multi-locus genotypes throughout localities to report the observed (Ho) and expected heterozygosity (He).